Comment on "Dynamical mean field solution of the Bose-Hubbard model" 



O 
(N 

>>■ 






I 

o 
o 



> 

o 

00 

oo 

o 
o 



^ 



K. Byczuk^ and D. VoUhardt^ 

'institute of Theoretical Physics, University of Warsaw, ul. Hoza 69, PL-00-681 Warszawa, Poland 
'"^Theoretical Physics III, Center for Electronic Correlations and Magnetism, 
Institute for Physics, University of Augsburg, D-86135 Augsburg, Germany 



In their preprint Anders et al. arXiv: 1004. 0510] propose a crucial modification of tfie Bosonic 
Dynamical Mean-Field Theory (B-DMFT) derived by us [Phys. Rev. B 77, 235106 (2008)]. Here 
we show that the modification consists of two steps which, in fact, cancel each other. Consequently 
their self-consistency equations are identical to ours. 



According to Anders et al. [l| they have modified 
the Bosonic Dynamical Mean-Field Theory (B-DMFT) 
y "to avoid double counting of kinetic energies" and 
thereby found stable solutions to the B-DMFT equations. 
Here we show that the alleged double counting is caused 
by an improper definition by Anders et al. [l[ of the con- 
densate wave function. In Ref. [l| the error introduced 
thereby is subsequently compensated by (mathematically 
uncontrolled) Bogoliubov shifts of the integration vari- 
ables. 

For simplicity we only consider the non-interacting 
case iU = 0) with nearest-neighbor hopping of the bosons 
on a lattice with finite coordination number Z since this 
is sufficient to illustrate our point. The generalization to 
finite U is straightforward. The B-DMFT maps the lat- 
tice problem onto a single site (say, i = 0) which is self- 
consistcntly coupled to two dynamical baths, one com- 
posed of normal and the other of condensed bosons [^ . 
For all coordination numbers 2 < Z < oo the exact par- 
tition function may be derived from a path integral with 
the local action Zexact = ^'■'''' /£'[&o&o] exp(-S'o[6o6o])- 
Here a superscript (0) refers to a lattice with a cavity, 
i.e., where the impurity site i = is removed, while a 
subscript refers to the impurity site itself. The local 
action for t/ = reads 



5o=y dTbl{T){dr - ^i)hQ{,T) 
K(f) / dTlb'^{Tl) + H.c. 
dn I dT2A(Tl - T2)b*{Tl)bo{T2), (1) 

where the condensate wave function 

= {b,)sm (2) 

and the hybridization function 

A(ti - T2) = Y.hntno{bn{ri)b*,{T2))sio, (3) 

are determined on the lattice with the cavity, and n ^ 

§■ t^j = Zt\sSi geometric factor; cf. Appendix A of Ref. 
for a derivation and notation. In particular we note 



that expectation values (. . . )^(o) only contain connected 
terms. 

In Ref. [l| a different approach is taken. Here the ex- 
pectation value of the boson is calculated on the impurity 
site 1 = 0, {bo)soj ^-nd this quantity is incorrectly identi- 
fied with the condensate wave function </>. By performing 
Bogoliubov shifts of the integration variables &o and Bq 
the corresponding error in (f> is compensated, whereby the 
geometric term k now acquires an additional, but spuri- 
ous term linear in the hybridization function, and thus 
takes the form k = Zt — J (ir2A(ri — T2). For this reason 
the authors of [l| refer to k as "the coupling between the 
impurity and the condensate" . 

It should be noted, however, that on a Bravais lattice 
with 2 < Z < 00 the two averages are not the same, i.e., 
(6j)g(0) 7^ {ba)so- This fact was already noted in Ref. [3| 
and was employed by them to solve the B-DMFT equa- 
tions numerically. In fact, as proved in the Appendix, 
the two local averages obey the following exact relation 
in Matsubara representation: 
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(4) 



where Gij are two-site correlation functions on the corre- 
sponding translationally invariant lattice. Using the cor- 
rect definition of 0, namely (f) = (6j)g(o) , and the identity 
I]ji *OjiGjio(w„) = -1 + ujnGooioJn) 3 it immediately 
follows that the spurious hybridization term in k intro- 
duced by the Bogoliubov shifts in Ref. jl| disappears 
again, whereby their self-consistency equations become 
identical to ours [2j. The claims by Anders et al. [l| 
concerning deficiencies of the original B-DMFT M are 
therefore invalid. In fact, their numerical results [l| are 
the solution of the B-DMFT equations previously derived 
by us in Ref. |3- 
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Appendix: Proof of eq. (|4]) 

The local correlation functions are determined via a 
functional derivative with respect to infinitesimally weak 
t/(l) symmetry breaking fields rji which are added to 
the action, i.e., {bj{T))g(o) = S\nZ^^^ /Srj*. The second 
derivative gives the two-site correlation functions 



Kiri)b%{T2))s^o^ 
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■{bji{ri))sw- (5) 
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In the linear regime of rji the solution of ([5]) reads 

(^ii(n))s(o) =J2j^{bjATi)b%{T2))smVj2- Employing the 
exact relation between two-site correlation functions de- 
fined on a lattice with a cavity, G\ { , and on a trans- 
lationally invariant lattice, Gj^^j^ , respectively, namely 



^3132 ~ '-'3132 
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we find 



(^ji)s(o) - / jGj^j^rij^ - Gj.qGqq Goj^Tfe) 



32 

^31 ~ '-'iioGgo 00, 



(6) 



where (i)j^ and ^o are local correlation functions at site ji 
and 0, respectively, on a translationally invariant lattice. 
In view of the translational invariance these quantities 
are the same, i.e., (pj-^ ~ c/jq, which proves Eq. ([4]). 
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